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ABSTRACT 



We present the data acquisition strategy and characterization procedures for 
the Canada- France Echptic Plane Survey (CFEPS), a sub-component of the 
Canada-Prance-Hawaii Telescope Legacy Survey. The survey began in early 2003 
and as of summer 2005 has covered 430 square degrees of sky within a few degrees 
of the ecliptic. Moving objects beyond the orbit of Uranus are detected to a mag- 
nitude limit of niR— 23 — 24 (depending on the image quahty). To track as large 
a sample as possible and avoid introducing foUowup bias, we have developed a 
multi-epoch observing strategy that is spread over several years. We present the 
evolution of the uncertainties in ephemeris position and orbital elements as the 
objects progress through the epochs. We then present a small 10-object sample 
that was tracked in this manner as part of a preliminary survey starting a year 
before the main CFEPS project. 

We describe the CFEPS survey simulator, to be released in 2006, which allows 
theoretical models of the Kuiper Belt to be compared with the survey discoveries 
since CFEPS has a well-documented pointing history with characterized detec- 
tion efficiencies as a function of magnitude and rate of motion on the sky. Using 
the pre-survey objects we illustrate the usage of the simulator in modeling the 
classical Kuiper Belt. 

^Based on observations obtained at the Canada-France-Hawaii Telescope (CFHT) which 
is operated by the National Research Council of Canada, the Institut National des Sciences 
de rUnivers of the Centre National de la Recherche Scientifique of France, and the Univer- 
sity of Hawaii, and on further observations obtained at the ESO Telescopes at La Silla and 
Paranal Observatories. In addition, based on data taken using ALFOSC, which is owned by 
the Instituto de Astrofisica de Andalucia (lAA) and operated at the Nordic Optical Tele- 
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1. Introduction 

Exploring the Kuiper belt allows us to study the process of planet formation and 
gain insight into the evolution of the outer Solar System. By studying the dynamical 
environment of these distant planetesimals, we can provide limits to planetesimal accretion 
and erosion models (Stern & Colwell 1997,b; Kenyon & Luu 1999). Clues to the outer 
Solar System's evolution are hidden in the relative populations of the various components 
of the trans-neptunian population (Malhotra 1995, 1996; Chiang & Jordan 2002; Levison 
& Morbidelli 2003; Gomes et al. 2005). A census of these dynamical populations allows 
us to identify potential source regions for the origin of short-period comets (Levison & 
Duncan 1997). However, an accurate picture of the Kuiper belt is difficult to determine in 
the presence of extreme observational biases related to flux and orbit determination. 

The biggest repository of Kuiper belt orbits is the Minor Planet Center^, whose lists 
include all reported TNOs with linked observational arc- lengths of two nights or more. 

scope under agreement between lAA and the NBIfAFG of the Astronomical Observatory 
of Copenhagen, observations collected at the Centro Astronmico Hispano Alemn (CAHA) 
at Calar Alto, operated jointly by the Max-Planck Institut fr Astronomic and the Instituto 
de Astrofsica de Andaluca (CSIC), and observations from the Hale Telescope, Palomar Ob- 
servatory, operated as part of a collaborative agreement between the California Institute of 
Technology, its divisions Caltech Optical Observatories and the Jet Propulsion Laboratory 
(operated for NASA), and Cornell University. Also, from observations taken as a visiting as- 
tronomer, Kitt Peak National Observatory, National Optical Astronomy Observatory, which 
is operated by the Association of Universities for Research in Astronomy, Inc. (AURA) 
under cooperative agreement with the National Science Foundation. 

^The Minor Planet Center database is online at http://cfa- 
www.harvard.edu/iau/mpc.html 
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While there are over 1000 Trans-Neptunian objects (TNOs) in the Minor Planet Center 
(MPC) database, over 300 of these objects are currently considered "lost" as their ephemeris 
uncertainty is now over 1 degree. The reasons why some objects are lost and other objects 
are not vary — foUowup observations have generally been done on an ad-hoc basis, so 
foUowup success depends on telescope scheduling, the brightness of the TNO, and how 
closely the initial orbit estimate matches the actual orbit. The orbit assumptions can be 
an important "invisible" bias, as these generally tend to place newly discovered TNOs on 
orbits similar to previously known TNOs. In addition, the foUowup observations are never 
systematically documented, so one does not know if a given object was lost because it was 
never followed (in which case the orbit listed in the database may be correct) or attempted 
and not found (in which case the database's orbit is certainly incorrect\) . The combination 
of these factors leads to possible biases in the orbit distributions in the MPC database. 
Out of the ~1000, only about 300 MPC TNOs have orbits well-determined enough for use 
in dynamical models, but even these lack information on the discovery/recovery conditions. 
Published individual surveys with characterized detection conditions, including some 
foUowup details, contain only tens of objects at most, providing a very small sample size. 

TNOs are loosely classified according to their orbits, into resonant, classical, scattered 
disk and extended scattered disk populations (see Luu & Jewitt (2002); Gladman et al. 
(2002); Gladman (2002). The extended scattered disk population was only relatively 
recently discovered (Gladman et al. 2002). Notably, if 2000 CR105 had not been tracked to 
a longer arc within its initial opposition outside its discovery survey, it would have appeared 
as a normal scattered disk object. Thus, the relative population of extended scattered disk 
objects cannot be determined without a knowledge of the tracking rate of discovered TNOs 
and a tracking program free from orbital biases. It is also important to have a large sample 
size, in order to be able to measure the population fraction of rare groups such as the 
extended scattered disk, Neptune Trojans (Chiang & Lithwick 2005), or TNOs in distant 



- 5 - 



mean-motion resonances; these classes appear to form only tiny fractions of a flux-limited 
trans-neptunian catalogue. 

The population densities of TNOs in each mean-motion resonance with Neptune, and 
the distribution of objects within these, help constrain models of Solar System evolution. 
For example, the relative fraction of plutinos (in the 3:2 resonance) versus "twotinos" (2:1 
resonance) might be used to determine the speed and extent of an outward migration of 
Neptune (Malhotra 1995; Ida et al. 2000; Chiang & Jordan 2002). In the 2:1 resonance, 
the speed of Neptune's migration and the initial conditions of the planetesimal disk affect 
the likelihood of capture into either leading or trailing libration islands (Murray-Clay & 
Chiang 2005). The longitudinal asymmetries caused by a resonant orbit dramatically 
alter the chances of finding resonant TNOs when pointing at particular ecliptic longitudes; 
the pointing history and search conditions of the survey are crucial for de-biasing from 
observed population fractions to the intrinsic ratios, even when objects are not discovered. 
This survey work has been designed to avoid these issues, and provide an exceedingly 
well-understood survey of the Kuiper Belt that comes with the tools needed to compare 
models with the survey's results. 

1.1. The CFEPS survey 

The Canada-Prance Echptic Plane Survey (CPEPS) is a part of the Canada-Prance- 
Hawaii Telescope Legacy Survey (see http://cfht.hawaii.edu/Science/CPHTLS/ for more 
details). The Legacy survey is divided into 3 components, with each component serving 
more than one science goal. As of July 31 2005, the Very Wide component of the Legacy 
Survey (LS-VW) has covered approximately 430 square degrees near the ecliptic, in multiple 
filters and with a widely-spaced exposure timing. The CPEPS Kuiper Belt study is the 
primary science driver for the LS-VW, and drives the exposure timing for this component 
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of the survey. Each field is acquired in multiple filters to serve other science goals (including 
galactic structure, stellar populations, and high-redshift quasar searches). However, the 
filter sequencing and exposure times are driven by the desire to detect TNOs and provide 
some portion of the tracking observations needed to determine their orbits for the CFEPS 
study. The goal of CFEPS is to obtain an unbiased census of Kuiper belt's various 
components. The survey will be well-characterized in terms of completeness and magnitude 
limits at each pointing, but also in the orbital tracking of each detected TNO. With the 
CFEPS sample, we intend to provide an orbital database suitable for dynamical modeling 
of the Kuiper Belt. 

In order to refine the CFEPS strategy, we conducted a small presurvey prior to the full 
onset of the project. This paper describes the CFEPS presurvey and how it affected the 
overall project design. The presurvey covered a 6.5-square degree patch of sky and resulted 
in the discovery of 13 TNOs; ten of these were within the expected fiux limit of CFEPS, 
and were followed in an intensive observational campaign. In addition, we describe the 
CFEPS simulator for comparing model Kuiper belts against our survey observations. 



2. Observations 

While this team has had significant prior experience in detecting TNOs (Gladman et al. 
2001; Petit et al. 2005; Allen et al. 2001), the presurvey observations detailed here were 
conducted specifically as a test of our design for the CFEPS observation strategy and a trial 
of our methods for obtaining and reducing the data. The queue programming, observation 
strategy (particularly optimum observation timing), and our ability to subsequently track 
the discovered objects was of particular importance. 
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2.1. The Observational Plan 

Given the CFEPS goal of producing a characterized TNO orbital database, the 
discovery observations must be characterized for limiting magnitude and completeness 
levels, and the recovery efficiency should be unbiased by orbital assumptions. In order to 
help attain this goal, CFEPS has some recovery observations built into the survey plan. 
These recovery observations are not targeted to specific ephemerides and thus are not 
affected by orbit assumptions. TNOs may still fail to be recovered due to moving in front 
of a bright star, into a gap in the camera, shearing out of the field coverage or simply 
being below the limiting magnitude of the observations, but not directly because of orbit 
assumptions. The generalized observational sequence is as follows: 

• Discovery consists of a triplet of images. Each pointing is imaged three times, with 
each image separated by 45 minutes to 1 hour. These images are searched for new 
TNOs and characterized for completeness and magnitude fimits, using methods 
described in Petit et al. (2004). 

• Nailing is a single image of the discovery pointing, taken within four nights of the 
discovery triplet (either before or after) . The close proximity of the discovery triplet 
allows this single image to be used to identify (aka "nail" ) the objects found in the 
discovery images. 

• Checkup consists of another triplet of images, taken within a few months (before or 
after) of the discovery images. This is primarily intended to improve the ephemeris of 
the TNO for one- year recovery. 

• One-year Recovery observations are similar to discovery - a triplet of images on one 
night, combined with a single image within the same dark run. These observations 
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add to the total observed arc, generating a well-determined ephemeris, but the arc 
length is still insufficient to allow precise determination of the orbital parameters. 

• Orbit Classification observations are again similar to discovery or one-year recovery. 
With these observations at least two years after discovery, the orbit is finally 
well-determined (usually to within 0.01 AU in semi-major axis), permitting the 
dynamical classification of the TNO. 

This was considered the minimum tracking strategy for CFEPS. The presurvey 
provided the opportunity to determine the optimum timing for each of these steps in 
order to minimize object losses and maximize observational efficiency. Many of the TNOs 
discovered in the CFEPS presurvey received much more extensive tracking observations, 
relatively densely spaced in time. These TNOs serve as test cases of when observations are 
maximally useful and possible within the CFEPS survey format. 

2.2. The Presurvey Observations 

The presurvey TNO discovery images were obtained at the Canada-France-Hawaii 
Telescope on August 5, 2002 using the CFH12K camera and the Mould R filter. These 
discovery images consisted of a set of 20 different field pointings, each pointing imaged three 
times during the night. The CFH12K camera has a 28' x 42' field of view, which allowed us 
to cover a contiguous 6.5 square degree patch of sky. The layout of these field pointings is 
shown in Fig. 1. One two-minute image at each of the 20 pointings was taken, and then the 
entire set of 20 pointings was imaged a second time, and then a third time. This produced 
a 1.5-hour time separation between each image in the triplet. Given typical rates of motion 
(all known TNOs move between 1.5 and 5 "/hr at opposition) this spacing is more than 
adequate to detect the moving TNOs. This was the natural result of cycling through the 
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sequence of 20 pointings, with each image having a 2 minute exposure time followed by a 
2 minute readout time. A 2 minute exposure time for these discovery images ensured the 
overheads associated with the CFH12K camera readout did not result in a less than 50% 
duty cycle. The presurvey observations provided valuable insight into the efficiencies of the 
telescope pointing system, guide star acquisition, camera overheads, and overall observing 
efficiency for the types of sequenced observations that are part of the CFEPS project. 

After searching the fields with our moving object detection pipeline (Petit et al. 2004), 
thirteen TNOs were discovered. They ranged in brightness from rriji = 22.3 — 23.9 (see 
Table 1 for a list), and heliocentric distances estimated to be between 40 to 60 AU. Because 
exposure times were two minutes and the discovery conditions were good - average seeing 
was 0.7" with very clear skies - the faintest presurvey TNOs are fainter than the average 
expected from the main CFEPS survey (which uses 80-sec exposures in median seeing 0.9"). 

The exposure time and filter was the same for the nailing image as for the discovery 
triplet. The nailing image for these twenty fields was not obtained until August 8 for some 
fields, or August 9 for others, due to bad weather at CFHT. We found that this four-night 
delay to be shghtly problematic, as the larger change in TNO position makes objects harder 
to detect in the nailing image as their on-sky error ellipse becomes much larger. Three of 
the TNOs were not visible on their nailing images due to problems like these. 

Some foUowup observations for all 13 TNOs (including those which did not get nailing 
measurements) were obtained at a variety of telescopes, utilizing significant resources to 
ensure as many objects as possible were tracked with as high a cadence as feasible, especially 
within the discovery opposition. Table 2 contains details of the observational campaign. All 
thirteen presurvey TNOs were re-observed one month after discovery, at either the Nordic 
Optical Telescope, the VLT, and/or at CFHT. These checkup observations consisted again 
of a triplet of images. At CFHT, the checkup images were of the same pointings as the 
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initial discovery images (as is the case for CFEPS), but the other telescopes were targeted 
to the predicted position of each object. Telescopes 4-m and up (such as the VLT and 
CFHT) were needed for some of these 1-month checkup observations. After this point the 
three TNOs fainter than itir — 23.7 were dropped from the foUowup effort; the CFEPS 
TNOs would seldom be this faint and these objects had been tracked for one month before 
being dropped due to magnitude considerations only. The remaining ten TNOs were then 
followed through winter 2004, at CFHT, the NOT, the VLT, KPNO 4-m, ESO 2.2-m, and 
Palomar 5-m. Each of these recovery observations generally consisted of a triplet followed 
by a naihng image on a following night. Except for CFHT in June 2003 and August 2004 
(where the images were part of the CFHTLS and thus non-targeted), each of the later 
recoveries was the result of a targeted observation. Because the dropping of the three 
faintest objects was not made on orbital grounds, this small sample is bias free. 

2.3. Optimal CFEPS Observations 

The experience gained with the presurvey was applied to optimize CFEPS 
observations, in cadence, field layout, and exposure time. The optimal timing of each 
astrometric observation is a balance between reducing the rate of false candidates and 
confusion with main-belt asteroids, uncertainty in the position of the TNO, and determining 
accurate orbital elements for each TNO as rapidly as possible. This balancing act is most 
crucial in the first year of observations, from discovery to one-year recovery. 

2.3.1. Observing cadence 

The discovery triplets are taken at opposition, as the sky motion then allows separation 
between TNOs and asteroids. Each image is part of a larger contiguous "block" of between 
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ls and 20 square degrees, similar to the presurvey field layout and as shown in Fig. 1. 
By grouping pointings into contiguous blocks, we enhance the likelihood of recovering 
the previously-detected TNOs without the need for an accurate ephemeris prediction. 
This minimizes the role of the assumptions required to create an orbit estimate based on 
short-track TNO observations, which induce potential bias in a tracking strategy based on 
pointed recovery. Checkup and recovery observations cover the same sky pointings as the 
discovery triplets; by placing the fields contiguously a smaller percentage of objects move 
out of the block of pointings. The compact layout of the block also ensures easy cychng 
through fields at the telescope, with a minimum of overhead. 

The nailing image (of the same field pointing) is taken within the next four days so 
that the positional uncertainty of the object on this single image is small enough to identify 
the TNO with a high degree of certainty. The nailing observation is not strictly necessary 
in order to find the TNO again two months later, as the positional error by that time may 
only grow to ~ 100". However, we find that the nailing observation improves the 1-year 
ephemeris prediction; with only the discovery triplet and a checkup triplet two months 
later the 1-year ephemeris error can reach one degree, at which point the object could be 
considered "lost". The utihty of the nailing observation is illustrated in Fig. 2; with the 
nailing observation the ephemeris uncertainty at one year is a factor of two lower. 

The checkup observations could in principle be conducted at any time in the next four 
months before the objects disappear behind the Sun. One might think that waiting as long 
as possible is best, and indeed Fig. 2 shows that with a 4-month checkup, 1-year semi-major 
axis uncertainties are a factor 3 smaller and 1-year ephemeris uncertainties are reduced by 
an order of magnitude when compared with checkup observations conducted at 2 months. 
Unfortunately, practical considerations intervene. 

This is due to the nature of the checkup observations. The CFEPS check-up images are 
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not targeted recoveries but rather re-imaging the initial survey pointings (due to secondary 
science goals for the LS-VW outside the scope of this project, but this also helps avoid 
orbital bias in tracking). Since these pointings 4 months after opposition are available only 
for a short time at the start of a night and are moving < l"/hr, scheduling a survey like 
CFEPS becomes risky and difficult. One observation per night at the beginning of several 
contiguous nights would be required and this is very difficult to acquire consistently and to 
analyze because of non-linear TNO motion over sequential nights. With CFEPS needing to 
acquire check-up observations on 15-20 fields per dark run, waiting until four months after 
opposition becomes impractical. Waiting until three months after opposition has similar 
problems and the TNOs are near the stationary points in their orbits. Therefore, due to 
practical considerations outside the orbital determination, the CFEPS checkup observations 
are conducted two months after opposition. As the ephemeris error ellipse is much smaller 
at two months compared to four, this is also a more conservative observing strategy. 

One year recovery observations consist of another triplet, followed by a single image 
on a nearby night, similar to the original discovery sequence. These images are taken near 
opposition, again to avoid asteroid confusion. A two-month deviation from this timing, 
such as would be possible within the queue, does not make much of a difference towards 
reducing uncertainty in semi-major axis or positional error, but would make detecting 
the TNOs much more challenging. The one-year recovery observation reduces future 
positional uncertainty to the order of an arcminute. However, with a two-month checkup, 
the semi-major axis uncertainty is still between 1 and 10 AUs after this one-year recovery 
(see Fig. 2). A further observation at or after two years is required to reduce the orbital 
parameter uncertainties another order of magnitude, to the levels necessary for dynamical 
modeling or resonance family determination. 

This CFEPS observation strategy is designed to minimize TNO losses, and particularly 
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to minimize any orbital bias in which objects are recovered. An insidious problem in TNO 
surveys has been recovery biases due to orbital assumptions. If the object is searched for 
on a later date, and not recovered, is it lost because it is below the limiting magnitude, 
in front of a star or in a chip gap, or because its orbit is actually quite different than 
predicted, and thus has moved out of the field of view? Because the CFEPS checkup and 
recovery observations are non-targeted, cover large patches of sky and occur at a rapid 
time cadence to avoid large ephemeris uncertainties, these objects should not be lost due 
to assumptions in the orbit. The other factors in losing objects — chip gaps, bright stars, 
limiting magnitude problems — have no correlation with the TNO orbit. Some objects 
in the survey will be lost out of the field of view, as the CFEPS checkup and recovery 
fields must observe the same pointings as the discovery fields, due to the Legacy nature 
of the CFHTLS-VW. The fraction of discovered objects which are not tracked will be 
characterized for each block. 

2.3.2. Exposure times and filters 

With Megacam, the camera used for CFEPS, the readout is much shorter than with 
CFH12K; 30 seconds instead of 2 minutes. As such, we are able to use shorter exposure 
times without dropping below a 50% duty cycle. While we were able to follow 10 out of 
13 prcsurvcy TNOs for 2 years, three of the faintest TNOs were lost after one month. 
This illustrates the flux bias present in all TNO foUowup. In particular, it illustrates the 
difficulty of tracking R > 23.7 TNOs in variable weather conditions, even with 4-meter and 
larger telescopes! To avoid this, we change to 70 second exposures for the discovery triplet 
in CFEPS, with a g' filter. This ensures that almost all TNOs discovered are brighter than 
R — 23.5 (assuming V — R — 0.5) and much more easily trackable at a range of telescopes, 
while still providing a decent duty cycle at CFHT. It also generates a discovery rate of 
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approximately one object per pointing. The nailing image is taken in the same filter, with a 
slightly lengthened exposure time (90s) to avoid losing TNOs below the limiting magnitude 
if the seeing varies. 

Grouping field pointings into blocks provides a convenient way to cycle through each 
triplet. With the lower overhead of Megacam, the typical time spacing between images is 
only 50 minutes, instead of 1.5 hours. With typical seeing conditions of 0.8", this is still 
adequate to detect the motion of TNOs as distant as 150 AU. 

Checkup is done in i' with a longer exposure time (180s) to detect the faintest TNOs, 
assuming a g'-i' color of approximately 0.8. To provide imaging of the same fields in 
different filters for other projects, the one year recovery is performed in r', with 110 second 
exposures in the triplet plus nailing image. The range of filters is necessitated by secondary 
science goals of the LS-VW, such as stellar populations and galactic structure studies. 

3. Data Reduction 

Our automated pipeline(Petit et al. 2004) was used to search each triplet of images for 
moving objects. The discovery triplets are the images which were searched and characterized 
using the pipeline, with variants being used to help identify moving objects in checkup and 
recovery images. 

3.1. The Automated Pipeline 

The presurvey discovery observations from CFH12K were trimmed and bias corrected, 
then flat-fielded using a median of all the discovery triplets combined with strong rejection 
criteria to remove stars and galaxies. Each individual flat-fielded image was then implanted 
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with artificial TNOs and searched for moving objects, using two different techniques. Two 
catalogues of all objects detected in each image were made, one produced using Source 
Extractor and the other via a wavelet technique. Each of these catalogues for the triplet 
were then compared, and stationary objects removed. Remaining objects were checked for 
movements compatible with TNOs in the field, and candidate objects which appeared in 
both sets of catalogues (Source Extractor and wavelet) were shown to an operator. Further 
details on the moving object pipeline process can be found in Petit et al. (2004). 

In the CFEPS presurvey, thirteen real TNOs were confirmed out of this set of 
candidates. Each discovery triplet was searched for objects moving between 1 and 6 "/hr at 
angles between ±23° with respect to the echptic plane. 

The same procedure will be used to find moving objects in each image triplet taken 
in the LS-VW. The CFEPS TNOs are detected in a much larger field of view (36 CCDs 
covering 1 square degree instead of 12 CCDs covering 0.32 square degrees), but the 
techniques have been well tested and scale up to the larger data set nicely. With the 
exposure times as above, given the TNO luminosity function measured by Gladman et al. 
(2001), we expect to find approximately one TNO per square degree. 



3.2. Detection Efficiencies 

The detection efficiency of the discovery triplets was calibrated as a function of 
magnitude and rate. After implanting artificial moving objects and running the automated 
software, the resulting hst of detected objects is compared to the list of implanted objects. 
It is obvious that the magnitude of a moving object will affect how frequently it is detected, 
as a lower signal to noise object is more likely to escape the cataloging process within the 
moving object software pipeline. 
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The rate of motion of an object also affects the detection efficiency. Objects wfiidi 
are moving very rapidly can travel distances of order the seeing disk (or, a more typical 
threshold of half the seeing disk) during an exposure, leading to "trailing loss" as the signal 
from the object is spread out over a larger number of pixels. Efficiency losses can also 
occur for objects which are moving very slowly, due to misclassification of moving objects 
as stationary, or rejection as TNOs due to perceived nonlinearity in their motion. This 
"misclassification loss" affects low signal-to-noise objects more severely than bright objects, 
due to poorer centroiding in each image. At higher motion rates, the poor centroiding effect 
becomes less noticeable, as the error drops to a smaller percentage of the measured motion. 

3.2.1. Presurvey efficiency 

We implanted a total of 42286 artificial objects with random motions across all ranges 
searched in the presurvey, and with a flat magnitude distribution ranging from 21.0 to 
25.0. The number of artificial TNOs placed in each chip was approximately 150. Although 
trailing losses are not an issue for the presurvey, as each exposure was only two minutes 
in length, we have attempted to measure our detection efficiency in both rate and motion. 
Fig. 4. Misclassification losses were also negligible, due to the combination of good seeing, 
long time gap between exposures in the discovery triplet, and searching to a minimum rate 
of l"/hr. 

As shown in Fig. 4, the efficiency is the percentage of objects detected in each 
magnitude/rate bin compared to the number of objects implanted in the bin. Even though 
we implanted over 40,000 artificial objects into this relatively small area (6.5 square 
degrees), the resulting efficiency function is still somewhat coarse. Given that the rate 
effects are negligible in the presurvey, we have simply fit a smooth double hyperbolic 
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tangent function to the magnitude efficiency falloff: 
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where A, Rc, Wi and W2 are equivalent to the maximum efficiency, rollover magnitude, 
and equivalent widths respectively (Petit et al. 2005). For the presurvey, A — 0.906, 
Rc = 24.20, Wi = 0.615, W2 = 0.309. This places the 50% efficiency level at tur = 23.9 
with a very steep falloff beyond this magnitude. Due to the rate cutoffs in the search, the 
efficiency is zero for rates < l"/hr or > 6"/hr. 



As the CFEPS is completed in blocks, each block will have an average efficiency 
function in rate and magnitude. Trailing losses are ignorable for the short exposure times 
of the CFEPS discovery images, so we will provide a similar smooth double hyperbolic 
tangent magnitude efficiency function for each block. The detection efficiency tables will 
also be made available online for CFEPS fields. 



The orbital information for each of the TNOs discovered in this survey is listed in 
Table 1. The orbits and uncertainties are calculated using the Bernstein & Khushalani 
(2000) software, version 2.0. 

We did not detect any confirmed plutinos in this sample. The presurvey field coverage 
is only 20 degrees away from Neptune, placing the plutinos near aphelion and at the 
faintest end of their magnitude range. Thus the expected fraction in the sample is lower 
and so this lack of plutinos is not entirely surprising. We also did not confirm detections 



3.2.2. CFEPS efficiency 



4. The Presurvey Discoveries 
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of any Scattered Disk objects (SDOs), although one of the three objects which was lost 
was potentially an SDO (K02O01, with a high eccentricity and inclination). This leaves 
our ten well-tracked detections with semi-major axes between 42.6 and 47.5 AU, and 
distances between 39.1 and 50.6 AU. One of our TNOs, K02O12, appears to be in the 2:1 
mean-motion resonance, while two others are potentially in the 7:4 resonance. 

In order to demonstrate our survey simulator, we will consider all of these objects to be 
"classical belt" TNOs. This avoids the problem of adding resonance angles into the orbital 
model, and with such a small sample, is unlikely to make a significant difference. And so, 
we now describe in the final section, the survey simulator. 

5. The CFEPS Survey Simulator 

One of the primary goals of CFEPS is to be able to compare the results of our 
observational campaign against theoretical models with high precision. Therefore much of 
the design of CFEPS has been focussed on the need to have sufficient information available 
so that a model of the Kuiper Belt can be 'passed through' a simulator of the CFEPS 
survey. These modeled detections can then be compared against the objects actually 
found. A statistical comparison of object properties can then be made against any available 
known parameters, such as magnitudes, distances, and orbital elements. In particular, 
knowledge of the pointing history of the telescope and the detection efficiency of each block 
of the survey is critical. If the magnitude and areal coverage is not uniform, different 
populations of Kuiper Belt objects are favored (for example, plutino detection is much 
poorer near conjunction with Neptune than away). An unbiased TNO orbital tracking 
program, as described above, is also critical. Otherwise dynamical classes could be missed 
or their relative populations changed (for example, if Extended Scattered Disk Objects were 
preferentially lost due to incorrect orbit assumptions). 
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The content of the CFEPS simulator will evolve as more and more blocks are added to 
the CFEPS survey; blocks will be added as their objects come to orbits with arcs 2 years or 
longer (the first release will thus occur in late 2006). In this section we present the general 
algorithm of the CFEPS simulator but apply it only to the small pre-survey sample. 

Because the orbital parameter space of the Kuiper Belt is vast, 'inverting' the observed 
objects to an intrinsic orbital distribution is problematic. Bottke et al. (2000) discuss 
this in the case of modeling the near-Earth object (NEO) detections from the Spacewatch 
telescope's NEO survey in order to determine the fractional contribution of several main 
belt input sources. The Kuiper Belt problem is somewhat more difficult due to the fact 
that we are directly modeling the source region itself which thus has many more parameters 
than the simple single number corresponding to a fractional weighting, as in Bottke et al. 
(2000). That is, we must model the orbital distribution interior to the classical belt itself, 
rather than just assign a single weight to the 'classical belt'. 

In addition, the importance of resonant interactions in the observed orbital parameter 
space requires that the constraints (coming from the resonant angles) between the various 
longitudes of the particles relative to the planets (mostly Neptune) be respected; TNOs in 
resonances come to perihelion (where the strong detection bias favors discovery) at certain 
ranges of ecliptic longitudes. In this paper, we will concentrate on modeling the classical 
belt, a much simpler problem. 

5.1. Simulator algorithm 

The CFEPS simulator begins with a model distribution of the Kuiper belt (or Kuiper 
belt region of interest) and applies a set of detection constraints related to the known 
history of the CFEPS survey (field location, magnitude depth, rate sensitivity) in order to 
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draw a user-determined number of model detections. These model detections can then be 
compared to the real detections to determine if the model can be rejected statistically. 

The Kuiper belt model consists of a file of one of several forms. The simplest consists 
a (typically large) number of (a, e, i) triplets, with a specification that their H absolute 
magnitudes should be drawn from a power law distribution of the form 

N{H <h)(X 10"^ Hmin <h< Hmax 

where i^min and Hmax are the minimum and maximum values and a is the slope of the 
power law. Published values for the slope a in the apparent magnitude range down to 
ruR ~ 24 range from 0.6 < a < 0.9 (Trujillo et al. 2001; Gladman et al. 2001; Elliot et al. 
2005; Bernstein et al. 2004). 

A second possible form for the model distribution allows the user to specify the H 
magnitude of each (a, e, i) triplet, although this will usually necessitate a much larger model 
file (to accurately represent the magnitude distribution). 

A third possible model distribution is one consisting of a six-element set (including 
the longitudes of node argument of pericenter lu;, and mean anomaly M) where the 
if-magnitudes are distributed by the same power law as above; this form is necessary 
at present if resonant angles are to be taken into account as the mean anomaly cannot 
be chosen independent of Vt and uj (note that Neptune's position needs to be taken into 
consideration). Lastly, the H magnitudes may be explicitly added to each of the 6-element 
sets instead (allowing different size distributions for different populations to be taken into 
account) . 

The CFEPS pointing history consists of a file with details corresponding to each block 
of the survey. Each block consists of a characterized region of search, parameterized by its 
width and height (in degrees), the RA and declination center of the block, the Julian date 
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of the pointing, the fiUing factor of the block (what fraction of the area inside the width and 
height specified is not in chip gaps for the Mosaic), the fraction of the discovered sample in 
that block which was tracked to >2-year arcs, and a pointer to the efficiency function for 
that block. 

The efficiency function for each block contains a definition of the magnitude efficiencies, 
either given as an analytic fit to the efficiency function (for example, Eqn. 3.2.1 for the 
presurvey) or as a lookup table. The rate cuts of the detection pipehne can also be specified. 
For example, the presurvey discovery data was searched only between retrograde rates of 
1 and 6 "/hour. The lower rate cut is important in accurately determining the ability of 
surveys to detect large distant objects in the scattered or extended scattered disks. 

The survey simulator then accepts a request for a certain number of model detections, 
and proceeds to iterate until this is successful or 10 biUion iterations have occurred. The 
simulator first determines the faintest H magnitude that could be seen in the faintest of the 
input blocks. The looping algorithm is: 

• Choose a random line from the model distribution. If necessary, determine a random 
H magnitude from the specified power law. 

• If the H magnitude is greater than the faintest possible, redraw. 

• If not defined, determine random values for fl, uj, and M, uniform on [— 7r,7r]; if these 
are not uniform an explicit distribution should be part of the input model. Calculate 
the apparent R-band magnitude of the object, given its distance and H. 

• Determine if the object is both in the field of view and within the rate cuts of the 
survey; if not, redraw. 

• Determine if the object is seen, using the magnitude efficiency function and a random 
number from 0-1 to simulate the detection probability; if not seen, redraw. 
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• Based on the 'tracking fraction', randomly decide if the object is tracked, but record 
non-tracked detections as well. 

• Record the orbital elements and H magnitude of the detection, along with apparent 
magnitude, distance, and sky coordinates. 

As one can imagine, when the sky coverage is small or the H distribution very steep, 
many random draws from the model distribution must be made to generate sufficient model 
detections. As the sky coverage and magnitude reach of the surveys grows, fewer trials 
are required. Depending on the input Kuiper belt models, the amount of time required 
to generate a "simulated survey" can vary widely. Generating 10,000 detections for the 
presurvey data requires about 1 hour of CPU time on a 3GHz-class PC for a simple 
case where objects are always bright enough to be seen in the presurvey field (but the 
random object does not necessarily fall in the presurvey field of view). A simple power law 
distribution for H under the same a, e and i conditions requires about 200 hours of CPU 
time. 



5.2. Survey Simulator Examples 

The survey simulator recreates the biases inherent in the survey, allowing users to 
thoroughly understand how the theoretical model translates to detected TNOs. We will 
now illustrate the use of the survey simulator and compare the presurvey objects to classical 
Kuiper belt models, as an example of the future use of the simulator rather than with the 
expectation of making strong conclusions with our 10-tracked-object sample size. 
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5.2.1. Illustrative Examples 

Let us consider a very simple model with each Kuiper belt object having a fixed 
semi-major axis a = 44 AU, and inclination i = 2.0° (as the presurvey fields were centered 
on an ecliptic latitude of 2°) . In order to illustrate the effects of survey limiting magnitude 
on detected TNO orbital parameters, we create a series of models where all objects have the 
same eccentricity, e =[0.01, 0.1, 0.2, 0.3], and a power-law absolute magnitude distribution, 
N oc 10°^, a =[0.6, 0.8, 1.0], 2.0 < H < 11.0. For these models, we have chosen to pick 
1,000,000 "detections" from the simulator in each run — far above what would actually 
be discovered in our small field, but useful to increase the signal to noise in the orbital 
parameters found by the simulated survey. 

Looking at the detected distance and mean anomaly distributions produced by the 
simulations of this toy model (Fig. 5), we can quantify some interesting effects. As the 
eccentricity rises, there is a stronger bias towards finding the simulated TNOs at perihelion. 
This is visible in the mean anomaly as an "aphelion hole" where dramatically fewer objects 
are detected near mean anomaly 180°. For a given eccentricity, the depth of the aphelion 
hole depends on the slope of the H distribution; thus determination of the eccentricity 
distribution would require a joint fit with the size distribution. The aphelion hole is 
also visible in the distance distribution as an increasing fraction of objects detected near 
perihelion. The strength of this perihelion bias is remarkable, and might be surprising 
if one considers that objects with eccentric orbits spend the majority of their time near 
aphehon. For our shallowest {a — 0.6) H-power law distribution of e = 0.3 TNOs, 49% of 
the simulated objects were discovered within 2 AU of perihelion, although they spend only 
12.7% of their orbital period within this distance. With a steeper (a = 1.0) power law, 
this fraction increases to 66% within 2 AU. This is caused by the steep falloff in TNO flux 
as distance d increases {oc d'^) combined with a cutoff in the detection efficiency at some 
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magnitude. Not all individual TNOs will exhibit this bias: for large TNOs which remain 
above the survey efficiency falloff, there is a larger peak in detection likelihood near aphelion 
due to the increased fraction of time spent at these distances during the orbit. However, 
for any distribution of TNOs where a significant fraction become fainter than the detection 
efficiency falloff at some point in their orbit, this strong perihelion bias overwhelms the 
aphelion peak. 

To explore this further, consider a second toy model with a = 44 AU, but with equal 
quarters of the model having the four eccentricities listed above, with a power law H 
distribution a — 0.8. Within each eccentricity grouping, there is still the same strong bias 
towards detecting objects near perihelion. Moverover, this bias towards finding objects at 
smaller distances is reflected in an overall bias towards detecting high-e objects. Although 
only 25% of the input toy model TNOs had e = 0.3, 44% of the detected objects were from 
the e = 0.3 group (with the rest of the 'detections' being 25% of e = 0.2, 17% of e = 0.1, 
and 15% of e = 0.01). Due to the power law H distribution, these detection ratios remained 
constant as the limiting magnitude of the simulated survey was altered, as long as that 
limiting magnitude did not surpass the limits of the input model magnitude distribution. 



5.2.2. Simulation of Presurvey 

We now compare our presurvey detections against three more complicated models, 
representing either empirical or theoretical models of the classical belt region. For an 
empirical view of the Kuiper belt, we will consider the Minor Planet Center database of 
classical TNOs as our model distribution. Due to the generally low eccentricities (e < 0.3), 
moderate distances, and easier observational tracking of objects in this region, it is tempting 
to attempt to use the group of objects in the MPC database as a representative sample 
of the true population. We therefore take the a, e and i of TNOs in the MPC database 
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with a between 40.5 and 47.5 AU and multi-opposition orbits as the orbital input to our 
simulator (304 objects, plotted in Fig. 6). As a theoretical view, we note that according 
to numerical simulations (Duncan et al. 1995), most of the classical belt with q >38 AU 
is stable for a > 40.5 AU. We thus build two other simple models by assuming that the 
classical belt is part of a disk that slowly declines in surface density (falling proportional to 
a~|) and was excited by some mechanism which left phase space uniformly filled in e; we 
populate each semimajor axis up to an e corresponding to g=38 AU, but weighted by 
since there is more phase space volume at higher e (see Fig. 6). This provides a Kuiper belt 
model with stable phase space evenly filled. Two versions of this model are considered, one 
with 40.5 < a < 47.5 AU (the "short stability" model) and one with 40.5 < a < 60.5 AU 
(the "distant stabihty" model), to test the presence of a cutoff near the 2:1 resonance at 
48.4 AU. 

While we would hope to be able to determine between the obvious difference of a 
model truncated at 47.5 AU and a similar model extended to 60 AU, we would also hke 
to distinguish between the short stability model and the MFC database of classical TNOs. 
The existence of an 'edge' to the Kuiper belt near 47.5 AU has been recognized in previous 
work (Dones 1997; Jewitt et al. 1998; Allen, et al. 2002; Trujillo et al. 2001; Trujillo & 
Brown 2001; Petit et al. 2005), but this has previously been based on statistical analysis 
of the distance distribution only. The existence and location of such an edge has important 
imphcations for models of the formation of the Kuiper belt (e.g. Hahn & Malhotra (1999); 
Gomes et al. (2004)). The difference between the MFC database and the theoretical 
model is more subtle (see Fig. 6); the MFC database shows a lack of low-e classical TNOs 
with g > 44 AU when compared to the stability model. It is possible that the lack of 
low-e detections in the MFC database is due to the decreased detectabihty of these high q 
objects, an effect which would be re-created in the survey simulator. The relative number 
of TNOs in this region can be used as a discriminator between models of the evolution of 
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the outer Solar System (e.g Gomes (2003); Levison & Morbidelli (2003)). 

We randomize the longitude of perihelion, the longitude of node, and the mean anomaly 
and pick 50000 objects from the simulator for each model within each group, setting the H 
power law distribution to be o; = [0.6, 0.8, 1.0]. Thus, we ran 9 simulations. 

The cumulative distance distributions from the simulations can be compared against 
the 10 presurvcy detections tracked for 2 years. We conducted a K-S test by calculating 
the maximum difference between the cumulative distance distribution of each simulation 
and the cumulative distance distribution of our survey. For each model, the likelihood of 
this difference was then cahbrated by drawing 1000 random samples of 10 objects from the 
50000 simulated survey discoveries, and calculating the differences between their cumulative 
distance distributions and the total simulation distribution. The results of our calibrated 
K-S test are shown in Table 3, and the a — 0.8 cumulative distributions of these models 
are shown in Fig. 7. Table 3 shows that our small sample does not rule out (P< 5%) any of 
these models based on the distance distribution alone, except for the 40.5 < a < 47.5 AU 
stability model with a — 1.0. 

The power of the survey simulator, when coupled with a quantitatively characterized 
survey, is in the ability to compare a wide range of model parameters against the true 
survey results. As such, we can look not just at the distance distribution, but also at 
eccentricity distributions, semi- major axis distributions, or the perihelion distribution. 
Table 3 contains K-S test results in these orbital parameters. The semi-major axis and 
eccentricity distributions are poor at constraining the two sets of models which are limited 
to semi-major axis between 40.5 < a < 47.5 AU, due to our limited sample size. However, 
the semi-major axis distribution expected from the distant stability model should have half 
of the detections with a > 50 AU, which is very different than the presurvey distribution 
where the largest semi-major axis was a — 47.7 AU. The eccentricity distribution is also 
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quite different, as the eccentricities of our presurvey objects are lower than the simulation 
results for a distant disk. This is due to the fact that we are mainly "discovering" TNOs 
with a > 50 AU and high-e when they are near perihelion in the simulation. It also is 
why the same model (40.5 < a < 60.5 AU) cannot be rejected based on the perihehon or 
distance distributions of our small sample, due to the strong perihelion bias explained in 
section 5.2.1. The perihelion distribution can be used to rule out the 40.5 < a < 47.5 AU 
stability model with a = 1.0 but this is influenced by our assumed g = 38 AU cutoff which 
does not account for TNOs in resonances which can have smaller perihelia nor the existence 
of scattered disk objects. 

We find that we can rule out a stable phase space model with 40.5 < a < 60.5 AU 
through the semi-major axis and eccentricity distributions of our presurvey objects. With 
a larger sample size, more stringent limits could be placed on the presence of a "cold" 
low-e disk outside the 2:1 resonance; as the perihelion distribution in Fig. 7 shows, almost 
all the a > 50 AU detections have periheha inside 50 AU - the low-e TNOs beyond the 
2:1 resonance will make up a very small fraction of the detected objects in a fiux-limited 
sample. 

The distance and perihelion distributions rule out a stable phase space model truncated 
at 47.5 AU with a steep TNO H distribution of o; = 1.0. However, we are unable to rule out 
a 40.5 < a < 47.5 AU stability model or an empirical MPC-based model with a shallower 
magnitude distribution. This is interesting, indicating the lack of low-e classical TNOs with 
g > 44 AU in the MFC may just be an observational artifact and the classical Kuiper belt 
could indeed uniformly fill all stable phase space. 

As a very simple test of what would be possible with a larger sample size, we simply 
cloned each of our presurvey objects and repeated the K-S probability testing. If a sample 
twice as large had the same d, g, a and e distributions, it would be possible rule out all 



-28- 



if-magnitude distributions of the MPC model based on the distance distribution, all the 
40.5 < a < 60.5 AU stability models based on the a and e distributions, and all except the 
a = 0.6, 40.5 < a < 47.5 AU stability model based on the perihelion distributions. 

6. The future of the CFEPS 

The CFEPS is being carried out using CFHT's Megacam, covering 1 square degree per 
pointing. One obvious benefit is the huge areal coverage possible in a short amount of time 
(fa 20 sq dcg in 3 hours to rUg = 23.5). This increases the discovery rate above previous 
survey levels, but it also increases the percentage of TNOs found in checkup and recovery, 
as the objects which shear out of the field of view become a smaller fraction of the total 
discoveries. In the first year of CFEPS results, we have found 65% of TNOs discovered 
within the survey can be followed to one-year arcs without using other telescopes. Using 
other facilities has allowed some blocks to obtain 100% recovery levels. 

We have settled on an observing strategy of discovery triplet, nailing, checkup at 
two months, followed by recovery at one year. The final measurement of semi-major axis 
and eccentricity will be obtained from a final recovery at either two or three years from 
initial discovery (depending on the whether the final observations for a given TNO arc 
are obtained at CFHT or another telescope). This appears to be an optimal balance of 
practical considerations and scientific gain. 

At present, over 100 TNOs have been discovered and tracked to one-year recovery 
within CFEPS, but these have still not been followed long enough to measure the orbital 
parameters accurately. With the presurvey results presented here, we are able to rule out a 
simple model of a uniformly filled stable phase space from 40.5 < a < 47.5 AU {q > 38 AU) 
if the magnitude distribution follows a = 1.0 and any uniformly filled stable phase space 
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models extending 40.5 < a < 60.5 AU with any iZ-magnitude distribution. Using the survey 
simulator, we are able to compare orbital distribution models with survey results, allowing 
detailed testing of models. By 2007, the CFEPS will be capable of detailed model testing, 
and we expect to be able to rule out the majority of models in the literature, allowing the 
community to refine the remaining (and newly proposed) models of the formation of the 
outer Solar System. 
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Fig. 1. — Field layout for presurvey (solid outlines) and CFEPS (dotted outlines) discovery 
images. The upper dashed diagonal line through the plot is the location of the ecliptic plane, 
the lower dashed diagonal is the location of the invariable plane. The locations of objects 
detected within the presurvey are marked on each field. The true aspect ratio is such that 
the CFEPS Megacam fields are square. 
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Fig. 2. — Time evolution of positional and semi-major axis uncertainty for two of our presur- 
vey classical-belt objects. Symbols are as described in Fig. 3. Panels (a) and (b) show the 
difference between two month and four month checkup. Panels (c) and (d) are similar, for a 
different object which did not receive a nailing observation. The uncertainties are generally 
low for these objects, as they are on nearly circular orbits. 
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Fig. 3. — Time evolution of positional and semi-major axis uncertainty, for extended scat- 
tered disk object 2000 CR105. On a log-normal scale, solid lines indicate the growth of posi- 
tional uncertainty with time, extended with a dotted line after further observations (at the 
vertical lines). The measured errors between the expected position and the actual position 
are indicated by the solid circles. On a normal-normal scale, the straight broken lines illus- 
trate the evolution of the semi-major axis uncertainty. Open squares show the a estimated 
uncertainty, as given by the Bernstein & Kushalani code. Solid triangles show the difference 
between the calculated semi-major axis at each observation and the final semi-major axis 
measurement. Note the much larger a errors compared with the previous examples (Fig. 2). 
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Fig. 4. — Detection efficiencies for tlie CFEPS Presurvey, split into bins of magnitude and 
TNO motion rate. Panel (a) is a viewpoint sliowcasing variation witli magnitude — note 
tlie steep drop-off at all rates near R = 24. This is consistent with previous TNO search 
efficiency calibrations. The view in panel (b) shows variation with rate (in pixels / hour) 
for each magnitude bin. Note the slight decrease in efficiency at low rates compared to the 
efficiency at a higher rate for similar magnitudes. 
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Fig. 5. — The mean anomaly and distance results of toy model simulator runs, with a range 
of eccentricities and H power-law distributions. Top: the mean anomaly distributions, with 
H distributions with a = 0.6 shown as solid lines and labelled with e. The dotted line 
{a — 1.0, e = 0.2) illustrates the degeneracy between the H distribution and the eccentricity 
distribution. Bottom: The distance distribution, a — 0.6 are shown as solid lines, a — 1.0 
are shown as dotted lines. The mean anomaly and distance distributions above show a bias 
towards detection at perihelion that strengthens with increasing eccentricity and steepening 
H slope. 
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Fig. 6. — Semi-major axes and eccentricities of the models used as input for the simulations 
of our presurvey discoveries. The solid circles are the objects from the MPC database used 
as an empirical model, with 40.5 < a < 47.5 AU and multi-opposition orbits. The small dots 
show the stable phase space models used, with 40.5 < a < 47.5 AU (short stabihty model) 
and 40.5 < a < 60.5 AU (distant stability model) subsets used as described in the text. 
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Fig. 7. — The cumulative distance, perihelion, semi-major axis, and eccentricity distributions 
for the presurvey objects (the stepped solid line) and the simulation results. The simulation 
objects are the results of running an MFC Classical Disk-based model with a. — 0.8 through 
the simulator (long-dashed line labeled "MFC" ), a stable phase space model with a — 0.8 
and 40.5 < a < 47.5 AU (dashed "Short Stability" line), and an extended stable phase space 
with OL = 0.8 and 40.5 < a < 60.5 AU (dotted "Distant Stabihty" line). 
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Tablc 1. Objects DisccA'crccl in tlic CL'EPS Prcsiirvc}'. 



Name Mag (R) (AU) a(AU) e i (°) Node(°)^ Arg of Peri. Time of Peri. 



K02O20 


22.3 


42.2 


42.648 


0.016 


1.271 


88.37 


291 


2466784 








± 0.009 


± 0.001 


± 0.001 


± 0.04 


± 5 


± 1272 


K02O32 


22.5 


47.4 


44.88 


0.074 


3.933 


339.814 


206 


2410815 








± 0.02 


± 0.002 


± 0.001 


±0.004 


± 2 


± 734 


K02O42 


22.6 


40.6 


43.63 


0.124 


5.450 


336.79 


289.24 


2437819 








± 0.02 


± 0.001 


± 0.001 


± 0.001 


± 0.30 


± 47 


K02P32 


23.1 


42.8 


42.65 


0.041 


1.570 


119.105 


112 


2426450 








± 0.03 


± 0.004 


± 0.001 


± 0.023 


± 1 


± 354 


K02P41 


23.2 


47.7 


44.34 


0.077 


2.510 


109.59 


28 


2502391 








± 0.01 


± 0.002 


± 0.001 


± 0.020 


± 6 


± 1818 


K02P12 


23.2 


50.6 


46.25 


0.144 


3.703 


328.776 


134.4 


2492063 








± 0.01 


± 0.001 


± 0.001 


± 0.000 


± 0.5 


± 161 


K02O03 


23.5 


39.09 


43.72 


0.191 


0.751 


182.71 


77.6 


2438640 








± 0.02 


± 0.001 


± 0.001 


± 0.04 


± 0.1 


± 12 


K02O12 


23.5 


47.2 


47.75 


0.221 


1.918 


329.086 


97.48 


2477337 








± 0.02 


± 0.001 


± 0.001 


± 0.001 


± 0.06 


± 33 


K02O43 


23.5 


42.7 


44.702 


0.059 


3.574 


336.730 


309 


2440637 








± 0.007 


± 0.001 


± 0.001 


± 0.002 


± 1 


± 278 


K02O40 


23.7 


43.0 


43.015 


0.030 


3.017 


110.81 


125.0 


2426965 








± 0.009 


± 0.003 


± 0.001 


± 0.01 


± 0.3 


± 95 
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Table 1 — Continued 



Name 


Mag (R) 


(AU) 


a(AU) 


e 


in 


Node(°)^ Ar 


g of Peri. 


Time of Peri. 


K02O41 


23.8 


42 


48 


0.2 


1.27 


69 


190 


2436197 






2" 


± 25 


± 0.6 


± 0.01 


± 3 


± 100 


± 20529 


K02O10 


23.8 


42 


46 


0.5 


3 


120 


317 


2469290 






7 


± 20 


± 1 


± 2 


± 23 


± 102 


± 12142 


K02O01 


23.9 


59 


40 


0.2 


1.0 


10 


130 


2494782 






± 5 


± 15 


± 0.5 


± 0.2 


± 12 


± 100 


± 22379 



^Distance from Earth. 

''These orbital parameters are as reported by Bernstein & Khushalani (2000). 

■^Uncertainty in D is only listed for objects with observed arc < 30 days. Other AD are < 
0.003 AU. 



Tal)lo 2. Objec-t Trac-kiiig 



Object 


CFHT^ 

8/02 


C.A. 

9/02 


NOT 

9/02 


VLT 

9/02 


CFHT 

9/02 


C.A. 

10/02 


Pal 

11/02 


CFHT 

11/02 


CFHT 

6/03 


ESO 

7/03 


VLT 

10/03 


CFHT 

8/04 


KP 

9/04 


Pal 

9/04 


K02O20 
K02O32 
K02O42 


5 

5,9 
5 


3,4,8 
8,9 
3,6,7 


3,4 
4 

3,4 


- 
- 


30 
30 
30 


3 

3,4 
3 


- 
- 


28,29,30 


23 
23 
24,25 


28,29 
28,29 
28,29 


- 
- 


12 


6,7 
6,7 


- 


K02P32 
K02P41 
K02P12 


5,9 
5,9 

5 


6,7 


4 
4 




30 
30 

30 






28,29,30 


24,25 
23 


29 
29 

28,29 




- 
12 


7 
11 


- 1 

to 

16 1 


K02O03 
K02O12 
K02O43 


5,8 
5,8 
5,9 




4 
4 


- 
4 


offW 
30 
30 


- 


5,6 








28,29 

29 
28,29 


12 
12 


11,12 




K02O40 
K02O41 
K02O10 

K02O01 


5,9 
5,9 
5,8 

5,8 






4 
4 
4 
4 


offS 
gap 
offW 
ofFW 












29 






16,17 



^Telescope abbreviations as follows: CFHT = CFHT 3.6-m + 12K camera, C.A. = Calar Alto 2.2-m, NOT = Nordic 
Optical Telescope 2.5-in + MOSCA, VLT = VLT + FORSl, CFHT (after 2003) = CFHT + Megacam, Pal = Palomar 5-in, 
ESO = ESO 2.2-m + WFI, KP = Kitt Peak 2.2-m and 4-m + Mosaic. 



I 

CO 

I 
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Table 3. Calibrated KS Test Results from comparison of simulated survey "detections" 
and presurvey detections. Sample size is 10 objects. 



Parameter 


H power law Index 
a = 0.6 a = 0.8 a = 1.0 


Distance 
Perihelion 
Semi-major Axis 
Eccentricity 


MPC model 40.5 <a< 47.5 AU 
18.2% 14.7 % 11.3% 
72.1% 81.1% 66.4% 
79.5% 71.6% 58.4% 
63.6% 66.9% 75.3% 


Distance 
Perihelion 
Semi- Major Axis 
Eccentricity 


Stability model, 40.5 <a< 47.5 AU 
45.6% 15.2% 2.4% 
14.2% 10.8% 1.5% 
51.8% 64.1% 62.4% 
55.5% 42.9% 13.8% 


Distance 
Perihelion 
Semi-major Axis 
Eccentricity 


Stability model, 40.5 <a< 60.5 AU 
61.6% 87.0% 40.2% 
18.5% 22.5% 36.2% 
0.0% 1.4% 0.0% 
1.5% 1.4% 3.0% 



